Dynamic Scaling in the Susceptibility of the Spin-| Kagome Lattice 
Antiferromagnet Herbertsmithite 
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The spin-i kagome lattice antiferromagnet herbertsmithite, ZnCu3(OH)eCl2, is a candidate ma- 
terial for a quantum spin liquid ground state. We show that the magnetic response of this material 
displays an unusual scaling relation in both the bulk ac susceptibility and the low energy dynamic 
susceptibility as measured by inelastic neutron scattering. The quantity xT a with a ~ 0.66 can be 
expressed as a universal function of H/T or ui/T. This scaling is discussed in relation to similar 
behavior seen in systems influenced by disorder or by the proximity to a quantum critical point. 

PACS numbers: 75.40.Gb, 75.50.Ee, 78.70.Nx 
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A continuing challenge in the field of frustrated mag- 
netism is the search for candidate materials which dis- 
play quantum disordered ground states in two dimen- 
sions. In recent years, a great deal of attention has 
been given to the spin- ^ nearest-neighbor Heisenberg an- 
tiferromagnet on the kagome lattice, consisting of corner 
sharing triangles. Given the high frustration of the lattice 
and the strength of quantum fluctuations arising from 
spin-i moments, this system is a very promising can- 
didate to display novel magnetic ground states, includ- 
ing the "resonating valence bond" (RVB) state proposed 
by Anderson A theoretical and numerical consensus 
has developed that the ground state of this system is 
not magnetically ordered [2]-l8|, although the exact ground 
state is still a matter of some debate. Experimental stud- 
ies of this system have long been hampered by a lack of 
suitable materials displaying this motif. 

The mineral herbertsmithite @, [ltj, ZnCu 3 (OH) 6 Cl 2 , 
is believed to be an excellent realization of a spin- 1 ; 
kagome lattice antiferromagnet. The material consists 
of kagome lattice planes of spin-i Cu 2+ ions. The su- 
perexchange interaction between nearest-neighbor spins 
leads to an antiferromagnetic coupling of J = 17±1 meV. 
Extensive measurements on powder samples of herbert- 
smithite have found no evidence of long range mag- 
netic order or spin freezing to temperatures of roughly 
50 mK 
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13]. Magnetic excitations are effectively gap- 
less, with a Curie-like susceptibility at low temperatures. 
The magnetic kagome planes are separated by layers of 
nonmagnetic Zn 2+ ions; however, it has been suggested 
that there could be some site disorder between the Cu 
and Zn ions [13. Ilij . This possible site disorder, with 
~ 5% of the magnetic Cu ions residing on out-of-plane 
sites, as well as t he p resence of a Dzyaloshinskii-Moriya 
(DM) interaction [16] , would likely influence the low en- 
ergy magnetic response. 

In this Letter we report a dynamic scaling analysis 



of the susceptibility of herbertsmithite as measured in 
both the bulk ac susceptibility and the low energy dy- 
namic susceptibility measured by inelastic neutron scat- 
tering. In particular, we find that the quantity \F a can 
be expressed as a universal function in which the en- 
ergy or field scale is set only by the temperature. This 
type of scaling behavior, when measured in quantum 
antiferromagnets 17 1 and heavy- fermion metals [181], has 
long been associated with proximity to a quantum criti- 
cal point (QCP). Power law signatures in the susceptibil- 
ity have also been associated with random systems such 
as Griffiths phase [3] or random singlet phase [2^] sys- 
tems. Such similarities could shed light on the relevant 
low energy interactions in herbertsmithite. 

Figure QJa) shows the ac magnetic susceptibility of a 
herbertsmithite powder sample as measured using a com- 
mercial ac magnetometer (Quantum Design). An oscil- 
lating field of 17 Oe, with a frequency of 100 Hz, was 
applied along with a range of dc fields up to hqH — 14 T. 
The data were corrected for the diamagnetic contribution 
by use of Pascal's constants. These results, for data sets 
with nonzero applied dc field, are plotted in Fig. [ljb) 
with x'T a (with a = 0.66) on the y axis and the unitless 
ratio hbH /ksT on the x axis. For this value of a, the 
data collapse quite well onto a single curve for a range of 
fiBH/ksT spanning well over two decades. Scaling plots 
with various exponent choices support a = 0.66 ± 0.02. 
This scaling remains roughly valid up to moderate tem- 
peratures, dependent upon the applied dc field. In the 
data taken with fioH — 0.5 T, shown in Fig. Q] the scal- 
ing fails for temperatures greater than roughly T — 35 K; 
under an applied field of [IqH = 5 T the scaling re- 
mains valid to about T = 55 K. The functional form 
of this collapse is qualitatively similar to the generalized 
critical Curie- Weiss function seen in the heavy-fermion 
compound CeCus.gAuo.ijUl, but with deviations demon- 
strating that such a simple response function is not quite 
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FIG. 1: (color online) (a) The in-phase component of the ac 
susceptibility, measured at 100 Hz with an oscillating field of 
17 Oe. (b) A scaled plot of the ac susceptibility data measured 
at nonzero applied field, plotted as x ac T a with a = 0.66 on 
the y axis and \ibH jksT on the x axis. Inset: A scaled plot 
of the dc magnetization, showing MT~ 34 vs fisH/hsT. 

adequate. It should also be pointed out that in herbert- 
smithite the entire bulk susceptibility obeys this scaling 
relation, while in CeCus.gAuo.i it is only the estimated lo- 
cal contribution, xl(T) = [xt^T 1 -%(T = 0) -1 ] -1 , that 
obeys scaling. A susceptibility of this form will imply a 
similar scaling in the bulk dc magnetization of the sam- 
ple, with MT a_1 expressible as a function of H/T. As 
a complementary measurement, such a scaling is shown 
in the inset to Fig. [ljb). The dc magnetization was 
measured up to /iqH = 14 T at temperatures ranging 
from T = 1.8 K to 10 K, and is plotted as MT" 34 vs 
H B H/k B T. 

The inelastic neutron scattering spectrum of herbert- 
smithite was measured on the time-of-flight Disk Chop- 
per Spectrometer (DCS) at the NIST Center for Neu- 
tron Research. A deuterated powder sample of mass 
7.5 g was measured using a dilution refrigerator with 
an incident neutron wavelength of 5 A. Measurements 
were taken at six different temperatures, with roughly 
logarithmic spacing, ranging from 77 mK to 42 K. The 
scattering data were integrated over a wide range of mo- 
mentum transfers, 0.5 < Q < 1.9 A -1 , to give a mea- 
sure of the local response. The momentum integrated 
dynamic scattering structure factor, S(ui), is shown in 
Fig-UJa). Similar to previous reports on the neutron scat- 
tering spectrum of herbertsmithite[ll[, the data show a 



broad inelastic spectrum with no discernable spin gap 
and only a weak temperature dependence for positive 
energy transfer scattering. The negative energy trans- 
fer scattering intensity is suppressed at low temperatures 
due to detailed balance. The imaginary part of the dy- 
namic susceptibility is related to the scattering struc- 
ture factor through the fluctuation-dissipation theorem, 
X"M = S(uj)(1 - e- Ru / kBT ). The dynamic suscepti- 
bility can then be determined in a manner similar to 
that used previously [ill. For the two lowest tempera- 
tures measured, detailed balance considerations will ef- 
fectively suppress scattering at negative energy transfer 
for values of \Tilj\ > 0.15 mcV. Thus these data sets are 
averaged together and treated as background. This back- 
ground is subtracted from the T = 42 K data, for which 
the detailed balance suppression is not pronounced below 
\hu\ = 2 mcV. From this, x"(cj; T = 42 K) is calculated 
for negative cj, and the values for positive uj are easily de- 
termined from the fact that x"( w ) is an odd function of 
u). The dynamic susceptibility at the other temperatures 
is calculated by determining the difference in scattering 
intensity relative to the T = 42 K data set. It is rea- 
sonably assumed that the elastic incoherent scattering 
and any other background scattering are effectively tem- 
perature independent. The calculated values of x"(w) 
at all measured temperatures are shown in Fig. EJb). 
The T = 42 K scattering data and x"( w ) were fit to 
smooth functions for use in calculating the susceptibility 
at other temperatures so that statistical errors would not 
be propagated throughout the data; the smooth function 
of x"( w i T = 42 K) used in the calculation is also shown 
in the figure. 

The resulting values for x"( w ) follow a similar scaling 
relation as the ac susceptibility, where the ratio huj/k B T 
replaces fi B H/k B T. In Fig.[3]we show x"MT a66 on the 
y axis and the unitless ratio Tioj jk B T on the x axis. The 
scaled data collapse fairly well onto a single curve over 
almost four decades of huj/k B T. Here we have used the 
same exponent a = 0.66 that was observed in the scal- 
ing of the ac susceptibility. However, the error bars on 
the data allow for a wider range of exponents (a = 0.55 
to 0.75) with reasonable scaling behavior. The collapse 
of the x"( w ) data is again reminiscent of the behav- 
ior observed in certain heavy-fermion metals, including 
the shape of the functional form of the scaling function. 
Let us assume that x"(u)T a oc F(u/T). The heavy- 
fermion metal CeCus.gAuo.i displays a scaling [ill I22I ] 
that could be fit to the functional form F(uj/T) = 
sm[cvian- 1 (w/T)]/[(o;/r) 2 + 1] Q / 2 . A fit to this func- 
tional form is shown as a dashed blue line in Fig. [3] This 
simple form does not fit the herbertsmithite data well for 
low values of u/T. Other heavy-fermion metals 23|, |24|, 
display a scaling relation that can be fit to the functional 
form F{lu/T) = (T/uj) a tanh(uj/l3T); this functional form 
is similar to that used to fit the dynamic susceptibility in 
Lai qfiSrn ruCuOdf^. This functional form fits our data 
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FIG. 2: (color online) (a) Neutron scattering structure fac- 
tor S(lo), measured using DCS integrated over wavevectors 
0.5 < Q < 1.9 A" 1 , (b) The local dynamic susceptibility 
x"((jj), determined as described in the text. Uncertainties 
where indicated in this article are statistical in origin and 
represent 1 standard deviation. 

much better, shown (with fit parameter f3 = 1.66) as the 
dark red line in Fig. |3l This function is somewhat un- 
usual, in that for low values of uj/T it is proportional 
to (uj/T) 1 - 01 rather than the expected w/T[17|; of course 
such a dependence might be recovered at still smaller val- 
ues of oj. For larger values oioj/T, this curve approaches 
a power law dependence with x"(w) oc w~ Q . This is con- 
sistent with the low temperature (T = 35 mK) behavior 
of the dynamic susceptibility of herbertsmithite reported 
earlier [1 1| . 

Other works on kagome lattice systems have shown 
evidence for similar behavior of the susceptibility. The 
dynamic susceptibility in the kagome bilayer compound 
SCGO has been shown to display power law behavior [26| 
and has been fit to a form|27[ identical to that shown as 
the dark red line in Fig.Hwith a = 0.4. Both SCGO and 
BSGZCOjiiJ demonstrate anomalous power law behav- 
ior in their bulk susceptibilities. Also, an early dynam- 
ical mean- field theory study of a kagome RVB state [2^] 
predicted such a scaling of the dynamic susceptibility. 
A recent paper 30] on herbertsmithite found that S(lj) 
was roughly independent of both temperature and en- 
ergy transfer for values of to greater than 2 meV. This 
simpler uj/T scaling is different from what we measure 
here in the low energy susceptibility. 

Similar scaling has been reported in other quantum 
antiferromagnets, many of which are believed to be close 
to a quantum phase transition [17j. The neutron seat- 
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FIG. 3: (color online) The quantity x"(u)T a with a = 0.66 
plotted against huj/kBT on a log-log scale. The data collapse 
onto a single curve. The lines are fits as described in the text. 

tering results on the spin-glass Lai.96Sro.o4Cu04 show 
a scaling that at small values of oj worked best with 
a = 0.41 ± 0.05[17| , while at higher energy transfers the 
data followed a pure uj/T scaling [2o]| with a = 0. Further 
comparisons can be made to neutron results on various 
heavy-fermion metals with doping levels that place them 
near a transition to a spin-glass or antiferromagnetically 



ordered state [22H24 |31| . In addition, a scaling of the 
ac susceptibility similar to that reported here was seen 
in CeCu5.gAuo.i[2l| and Ce(Run sRhn .O^Si?}^. Recent 
exact diagonalization work[33j has suggested that the 
ground state of the spin-i kagome lattice antiferromag- 
net with a Dzyaloshinskii-Moriya interaction will be a 
quantum disordered state for the Heisenberg Hamilto- 
nian, but a Neel ordered state when the component of 
the Dzyaloshinskii-Moriya vector perpendicular to the 
kagome lattice plane exceeds « J/10. The presence of 
a nearby QCP is also possible in models without a DM 
interaction [34|. Futhermore, several of the theoretically 
proposed ground states for the spin-i kagome lattice 
antiferromagnet [35l |36| are critical or algebraic spin liq- 
uid states. These proposed quantum ground states would 
possibly display excitations that are similar to fluctua- 
tions near a QCP. Thus, the observed low energy scaling 
behavior in herbertsmithite might signify quantum crit- 
ical behavior [37] or a critical spin liquid ground state. 

In many doped heavy-fermion metals, the observed 
non-Fermi liquid behavior is likely related to disor- 
der. In herbertsmithite, the low temperature suscepti- 
bility roughly resembles a Curie tail, and it has been 
suggested [Til [l5| that this is attributable to S = 1/2 im- 
purities (consisting of « 5% of all magnetic ions) with 
weak couplings to the rest of the system. We find that 
the scaling behavior seen in herbertsmithite does have 
features in common with the disordered heavy-fermion 
metals, such as x"( u ) proportional to (w/T) 1_Q at low 
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values of lj/T rather than linear vaoj/T. The divergence 
of the low temperature susceptibility may also be indica- 
tive of a random magnetic system, such as a Griffiths 
phase HU or random singlet phase [20l]. A collection 
of impurity spins subject to a broad distribution of cou- 
plings, P(Jimp) may result in a power law susceptibility 
at sufficiently low temperatures [39(; however, the scaling 
presented here describes the entire measured susceptibil- 
ity rather than the response of a small impurity frac- 
tion. In this scenario, the scaling of the ac susceptibility, 
with x'T" ' 66 oc F(H/T), would be useful in determining 
the distribution of couplings experienced by the impurity 
spins in herbertsmithite. The disordered heavy-fermion 
metal Ce(Ruq^Rh . 5)2812 displays a scaling of the ac 
susceptibility |32j that is remarkably similar to that in 
herbertsmithite, except with a considerably larger effec- 
tive moment [a much smaller field will suppress the low 
temperature susceptibility of Ce(Ruo.5Rho. 5)2812]. That 
scaling was attributed to a broad distribution of coupling 
strengths that likewise diverges at low coupling, quite 
likely with a power law distribution [ioj]: P(Ji mp ) oc J^ p - 
In terms of our data, this would imply a distribution of 
impurity couplings which extend to several meV, which 
is surprisingly large for the assumed out-of-plane impu- 
rity ions. A further prediction in this scenario is that the 
distribution of local susceptibilities would diverge at low 
temperatures [l9| such that the width of the Knight shift, 
8 K, as measured by NMR would go as oc T~ A with 
A ~ 0.17. The observed NMR signal certainly broadens 
at low temperatures 41], and it would be most interest- 
ing to see if it follows this specific power law. Ther- 
mal transport measurements would be important to help 
differentiate between scenarios where the scale-invariant 
spin excitations are localized near impurities or extended 
(as in the aforementioned criticality scenarios). 

In conclusion, we have shown that the low energy dy- 
namic susceptibility of the spin-^ kagome lattice anti- 
ferromagnet herbertsmithite displays an unusual scaling 
relation such that %T a with a ~ 0.66 depends only on 
the thermal energy scale fc^T over a wide range of tem- 
perature, energy, and applied magnetic field. This be- 
havior is remarkably similar to the data seen in certain 
quantum antiferromagnets and heavy-fermion metals as 
a signature of proximity to a quantum critical point. In 
addition to scenarios based on impurities, the results may 
indicate that the spin-i kagome lattice antiferromagnet 
is near a QCP, or that the ground state of herbertsmithite 
may behave like a critical spin liquid. 
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